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Abstract 

We study the Tevatron signatures of promptly-decaying slepton co-NLSPs in the 
context of General Gauge Mediation (GGM). The signatures consist of trileptons plus 
Fj T and same-sign dileptons plus J£ T . Focusing first on electroweak production, where 
the Tevatron has an advantage over the early LHC, we establish four simple bench- 
mark scenarios within the parameter space of GGM which qualitatively capture all the 
relevant phenomenology. We derive limits on these benchmarks from existing searches, 
estimate the discovery potential with 10 fb _1 , and discuss ways in which these searches 
can be optimized for slepton co-NLSPs. We also analyze the Tevatron constraints on a 
scenario with light gluinos that could be discovered at the early LHC. Overall, we find 
that the Tevatron still has excellent reach for the discovery of SUSY in multilepton 
final states. Finally, we comment on the possible interpretation of a mild "excess" in 
the CDF same-sign dilepton search in terms of slepton co-NLSPs. 
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1 Introduction 

1.1 Motivation 

Signatures involving several high leptons and missing energy are generally considered to be 
among the most promising channels in which to search for new physics at a hadron collider. 
Experimentally, these signatures are very clean. Leptons (by which we mean electrons and 
muons, following standard experimental parlance) are straightforward to identify and difficult 
to fake. Moreover, the production rates for processes leading to several leptons plus missing 
energy are naturally small in the Standard Model, so backgrounds for new physics are low 
in these final states. 

On the theoretical side, these signatures are also very well motivated. New physics 
scenarios commonly give rise to multiple leptons and missing energy, especially in supersym- 
metric models. (For an overview of common supersymmetric models and their phenomenol- 
ogy, see [T].) Perhaps the best known example occurs within the "minimal supergravity" 
(mSUGRA) framework, in which the vast number of soft parameters of the MSSM are sim- 
plified down to just four continuous parameters. In this model, the lightest supersymmetric 
particle (LSP) is a bino-like neutralino, and it is absolutely stable. Pair producing heavier 
neutralinos and charginos and decaying them down to the LSP leads to final states with 
multiple leptons and missing energy. 

In this paper we study a scenario that is at least as well motivated theoretically, although 
perhaps less known to experimentalists, in which multilepton plus E T signatures also au- 
tomatically arise. This is gauge mediated supersymmetry breaking (GMSB) with slepton 
co-NLSPs. (For a review of the phenomenology of gauge mediation and original references, 
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Figure 1: An example slepton co-NLSP process leading to a trilepton final state. Here, a 
charged and neutral wino are produced, and they decay through a right-handed smuon and 
stau, to produce a final state with fjT h + t~ + J£ T . 

see 0.) 

Recall that in gauge mediation, the LSP is a nearly-massless gravitino. The lightest 
MSSM sparticle is then the next-to-LSP (NLSP), and it always decays in a universal way to 
the gravitino plus its SM superpartner. Since this decay rate is heavily suppressed by the 
SUSY-breaking scale, all heavier sparticles decay first down to the NLSP before decaying to 
the gravitino. Therefore, the nature of the NLSP is the most important aspect of the GMSB 
spectrum for collider signatures. If the three flavors of right-handed sleptons all dominantly 
decay to gravitinos, and are thus co-NLSPs, then all MSSM events will contain at least two 
high pt leptons or taus, plus missing energy from the gravitinos. Additional leptons can arise 
from cascade decays of heavier SUSY states, so that it is common to have three or more 
leptons plus taus in the final state (for an example see figure [I]). We focus on the situation 
where the sleptons decay promptly, which means gravitino masses in the sub-keV range. It 
is also interesting to consider delayed decays, but these signatures are entirely different and 
will be considered elsewhere. 

The importance of multilepton plus Fj T signatures for the slepton (co-)NLSP scenario was 
first pointed out in a number of the original phenomenological studies of gauge mediation, 
including Refs. [31 H El El El [9] , and by the Tevatron Run II SUSY working group [10]. In 
this paper, we revisit these signatures. This is motivated by two reasons. 

First, the previous phenomenological work was focused on Minimal Gauge Mediation 
(MGM), a simplified model of GMSB in which the vast complexity of the MSSM is again 
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reduced down to a handful of parameters. Recently in [11] [12] , a general framework of gauge 
mediation was formulated, and the full parameter space was shown to be much larger than 
that of MGM. While MGM enforces several large hierarchies among the soft masses (for 
example the left-handed sleptons are about twice as heavy as the right handed sleptons and 
the colored superparticles are significantly heavier than the uncolored ones), General Gauge 
Mediation (GGM) requires no such hierarchies, allowing new combinations of states to play 
important roles in collider signals. Thus it is interesting to revisit the collider signatures of 
GMSB within the full parameter space of GGM, in a model-independent way. 

Second, we now have the opportunity to confront the GMSB signals with the large existing 
dataset collected by the Tevatron. Unfortunately, there was not much experimental follow- 
up to the early phenomenological work on GMSB with slepton co-NLSPs. While these 
models were searched for systematically at LEP (we review the limits below), Tevatron 
searches optimized for GMSB have not been conducted in the multilepton plus E T final 
states. However, related searches motivated by other scenarios (such as mSUGRA) have been 
conducted for same-sign dileptons p31 E] an d trileptons [151 EE EE] • Thus it is interesting 
to examine the constraints these searches place on GMSB with slepton co-NLSPs, and to 
understand how to optimize their discovery potential. 

We emphasize that the Tevatron remains very relevant despite the advent of the LHC. 
With 10 fb _1 , the Tevatron can produce more uncolored SUSY states than the LHC with 
< 1 ftT 1 at 7 TeV, while many of the backgrounds are larger at the LHC. Thus it will remain 
the best place to search for the production of light uncolored states in the near term [18J. On 
the other hand, the LHC will very soon have the advantage for producing new light colored 
states. Therefore the two programs will remain complementary for the next several years. 

1.2 Summary of the main results 

Here we attempt to give a somewhat self-contained summary of our approach and our results, 
for the benefit of the more casual reader. Our philosophy and approach, motivated by GGM, 
is very much a continuation of the program begun in [T9l [20] , where neutralino NLSPs were 
thoroughly studied in a general, model-independent way. (For other related references on 
collider signatures of GGM, see (2U E21 1231 M> ESI EE].) 



4 



Even after restricting the NLSPs to sleptons, there are many possible spectra within 
GGM, and the parameter space is too large to fully explore. In order to reduce the complex- 
ity, our strategy is to choose spectra with as few light states as possible which play a role in 
the signal of interest, with other states decoupled. To that end, we focus on production of 
electroweak SUSY states - neutralinos, charginos and sleptons - and assume that the colored 
sparticles are too heavy to be produced. This is a natural assumption, given the relative 
strengths and weaknesses of the Tevatron vs. LHC described above. 

With the remaining electroweak SUSY states, we identify (in section [3]) four simple 
spectral types, classified primarily by the type of sparticle being produced. Specializing 
to gauge eigenstates, these correspond to: wino production, higgsino production, and left- 
handed slepton production with and without an intermediate bino. For each spectral type, 
we identify a two dimensional parameter space, spanned by the production mass and the 
NLSP mass, which together control the production cross-section and signal kinematics. We 
choose parameters such that electrons, muons, and taus are produced democratically. 

These simplified parameter spaces are intended to serve as benchmark scenarios which 
contain within them all the relevant Tevatron phenomenology of GMSB with slepton co- 
NLSPs. They are model-independent, in that they are defined using weak-scale soft pa- 
rameters, and they are not subject to any arbitrary constraints or relations that exist in 
e.g. mSUGRA or MGM. Nevertheless we can be assured that our benchmark scenarios do 
correspond to physical models, since they conform to the parameter space of GGM [12]. 
Using these benchmark scenarios, we can investigate the phenomenological possibilities of 
slepton co-NLSPs. They can be used to quantify the strengths and weaknesses of existing 
experimental searches, and they can guide in the planning of future searches. 

To illustrate this, we use simulations to estimate the current limits and future (i.e. 
10 fb" 1 ) discovery reach for each parameter space from the existing Tevatron searches in 



the trileptons+ E T and SS dileptons+ E T final states. ^ Our main results are summarized 
in figs. [7] and [8] of section |5j Apart from regions of parameter space where decays are kine- 
matically squeezed, our results roughly translate to limits on the mass of the sparticle being 

(i)We use Pythia6.4 [57] to generate events, PGS4 [28 for detector simulation, Suspect [5S] to calculate 
supersymmetric spectra, and a combination of Pythia and private code to calculate decay tables. For more 
details on our simulations, see appendix [X] 
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produced: 

™>w £ 190 GeV, mfj > 190 GeV, m ih > 155 - 170 GeV (current limits) 
m w > 220 GeV, m A > 250 GeV, m ih > 230 - 240 GeV (10 fb" 1 reach) 

We find that the trilepton searches set the best current limits, owing to their larger dataset 
(3.2 fbr 1 vs. 1 ffcr 1 ). However, we find that the CDF SS dilepton search has the best 
projected reach, owing to its more inclusive nature. 

Also in section [5j we show that the Tevatron still has 3a discovery potential for a wide 
range of superpartner masses consistent with the current limits, and 5cr discovery potential 
for a narrower range of allowed masses, especially if the searches are optimized for the 
kinematics of the slepton co-NLSP signal. 

In section [6j we go beyond our four benchmark scenarios to discuss various promising 
directions for future research. First, we consider the degradation in experimental sensitivity 
for models that preferentially produce taus instead of electrons and muons, such as models 
with stau NLSP. Overall, we find that the limit and reach for such models is currently quite 
poor, given the higher backgrounds and lower efficiencies. The trilepton searches do end up 
being more sensitive than the SS dilepton searches, because they include one-prong hadronic 
r decays in their list of final states. Second, as preliminary to studying colored production 
at the LHC, we present the Tevatron limits and reach for a slepton co-NLSP scenario with 
gluino production. In a forthcoming paper [30], we will study the phenomenology of this 
and other GGM scenarios in much more detail. 

Section [7] contains our concluding discussion. Here, we suggest ways to optimize the 
Tevatron multilepton searches for GMSB, guided by our results in the previous sections. 

• These searches have been biased in various ways to accommodate the kinematics of 
the mSUGRA signal. As we will see, GMSB with slepton co-NLSP produces leptons 
with larger p^'s and larger missing energy than mSUGRA (a comparison appears in 
appendix [C]) . Thus the cuts on these quantities can be easily tightened without losing 
much signal acceptance. This could also help improve the sensitivity in channels that 
include hadronic r decays (e.g. the single track category of the trileptons search). 

• In GMSB it is quite common to have more than three leptons (up to 6-8 depending on 
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the model!) in an event, coming from cascade decays of heavier states [U [7J [8]. Thus 
it is important to be as inclusive as possible and not veto e.g. on additional leptons in 
the event. 

• Similarly, one should avoid vetoing on leptons forming Z bosons if at all possible, 
since with so many leptons in an event, it is quite common that two will accidentally 
reconstruct a Z. 

• Finally, we point out that there are two major gaps in Tevatron coverage of slepton 
(co-)NLSPs. First, there are no existing Tevatron searches for > 4 leptons plus missing 
energy, which is known to be a promising channel for discovery of GMSB [8j. Second, 
searches for SS dileptons which include hadronic tau decays would greatly improve the 
coverage of tau-rich scenarios, despite the higher backgrounds. 

Finally, we explore (also in section [7]) the intriguing possibility that a hint of an excess in 
the 1 fb _1 CDF same-sign dilepton search [TJ] is due to slepton co-NLSPs. The slight "excess" 
(not statistically significant) occurred at high Fj t and also at high lepton p?. Specifically, 
CDF saw four events with Fj T > 80 GeV and four events with lepton > 80 GeV where they 
expected ~ 1. While these events most likely came from insufficiently understood background 
tails, it is entertaining to consider their possible origin in our slepton co-NLSP scenarios. As 
discussed above, events with high J£ T and high lepton are very characteristic of slepton 
co-NLSPs. Indeed, we find that if these excesses are real, slepton co-NLSPs can easily 
explain them consistent with current experimental limits, especially the slepton-production 
benchmark scenarios where the existing limits are much weaker. Moreover, if these excesses 
are real and are explained by slepton co-NLSPs, they will show up spectacularly in 10 fb _1 
as a 5cr discovery. 

Our appendices include technical details and special topics. Appendix [X] contains details 
of our simulations. In appendix |Bj we apply our general results to the special case of MGM 
with slepton co-NLSPs. This is contained within just one of our benchmark scenarios (slepton 
production with intermediate bino), and we derive the corresponding limits on the MGM 
parameter space. Finally, in appendix [C] we compare the multilepton signals of GMSB to 
those of mSUGRA. 



7 



2 Slepton co-NLSPs in GGM 



2.1 GGM Parameter Space and Production Modes 

We begin by reviewing the parameter space of GGM (TTJ [TJ]. In GGM, the spectrum is 
determined by the five sfermion mass-squareds, m|~, m|, m|, rrv~ L , m\ R , which are subject 



to two sum rules^ and three gaugino masses, Mi, M2, M3, not subject to any relations. 
The sfermion mass-squareds are flavor degenerate at the messenger scale, but the first two 
generations become split from the third via RG evolution. (For a nice summary of the 
MSSM RGEs, see pQ.) The Higgs soft masses rn 2 Hu and m 2 Hd are sensitive to corrections from 
additional sector(s) that e.g. address the \i problem (see [H] |3T] for discussions of this in 
the context of GGM). So we treat these as free parameters, and trade them in practice for 
fi and tan /3 defined at the weak scale. Finally, the A-terms are zero at the messenger scale 
(and hence small at the weak scale), and the gravitino mass is <C m Bm t. 

As discussed in the introduction, since we are focusing on the Tevatron in this paper, 
we will assume electroweak sparticle production only, with colored sparticles decoupled, 



m coiored <L 500 GeV, ( 3 ) So the relevant parameters are: 



Mi, M 2 , /x, tan/3, mg L , m gfl , m ?L , m fR (2.1) 

Keeping in the spirit of our model-independent philosophy, we choose to specify these pa- 
rameters at the weak scale and remain agnostic as to their UV (messenger scale) origin. 

Potential SUSY production modes involve the neutralinos and charginos, or the left- 
handed sleptons (both charged sleptons and sneutrinos). We specialize for simplicity to 
gauge eigenstate -inos. Plots of wino, higgsino, and slepton production cross-sections are 
shown in fig. [2j Higgsinos and winos have similar cross-sections, with a ~ 100 fb for ~ 
~ 200 GeV. The total slepton production cross-section is a factor of 5-10 lower at the 

( 2 ) The GGM sum rules (see Ref. [TT]) do not play an important role here. This is because, for simplicity, 
we consider spectra with a small number of light superpartners that produce the signal of interest, with the 
other superparticles decoupled, m > TeV. For the examples we consider, the decoupled superparticle masses 

can always be chosen such that the sum rules are satisfied. 

( 3 ) Near the end of the paper, we will also briefly consider the Tevatron limits and reach for colored sparticle 

production, with an eye towards discovery at the LHC. 
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Figure 2: Tevatron production cross-sections (LO) for winos, higgsinos, and left-handed 
sleptons. For all plots, we take tan/3 = 2. For left-handed slepton production, we take 
flavor symmetric soft terms at the weak scale, and we fix /i = 300 and m TR = 100 GeV. 
These parameters enter the stau mixing. For wino production, we fix /i = M\ = 1 TeV. For 
higgsino production, we fix Mi = M 2 = 1 TeV. 



same mass, e.g. it is ~ 10 fb for 200 GeV sleptons. Not shown here is the cross-section for 
right-handed sleptons, which is further suppressed. As we will discuss in section 



2.3 



'R*-R 



production leads solely to OS dileptons plus Fj t , and the large SM backgrounds for this final 
state, combined with the smaller SUSY cross section, make this an unsuitable search channel 
at the Tevatron. 



2.2 Slepton co-NLSP Definition 

Now that we have defined the parameter space, we are ready to discuss what it means for 
a spectrum to have slepton co-NLSPs. We specialize to spectra where the right-handed 
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sleptons are at the bottom of the spectrum/ 4 ) Examples of such spectra are shown in fig. 



[3j these will be discussed in more detail in section |3j Since the right-handed slepton masses 
originate from a flavor universal boundary condition at the messenger scale, the right-handed 
selectron and smuon are always very nearly degenerate. The lightest stau, fi, is always the 
lightest state, because of RG running and left-right stau mixing. 

Whereas f\ will always decay to a r and gravitino, the other sleptons have two possi- 
ble decay channels [I]. For smaller Smj = rrie R — m^, the first two generation sleptons 
dominantly decay to a lepton and gravitino, and this defines the slepton co-NLSP regime. 
In a collider, every SUSY event contains at least two hard leptons or taus (depending on 
whether Ir or f\ is produced), and significant missing energy, carried by the gravitinos. For 
larger 5m z ~ , the first two generations dominantly decay three-body, to f\ through an off-shell 
neutralino, 1^ — > ^N* — > / ± r =F ff t . Every SUSY cascade ends in fi, and this defines the stau 
NLSP regime. 

Both the slepton co-NLSP and the stau NLSP scenarios are "generic" in models of gauge 
mediation, i.e. no special fine tuning of UV parameters is required to realize either scenario. 
In general, the slepton co-NLSP regime applies for Smj R < 5 — 10 GeV, while the stau 
NLSP regime holds for larger mass splittings |32j. The mass splitting in turn is controlled 
primarily by tan /3 and the left-right mass splitting, as this sets the amount of left-right 
mixing for the third generation sleptons. For example, in MGM tan /3 < 10 — 15 corresponds 
to dmf < 5 — 10 GeV [5J. Since we will focus on the slepton co-NLSP regime, we will 
henceforth take 

m- eR = m fR m- eL = m fL (2.2) 
This simplifies the parameter space, without much loss of generality. 



2.3 Prequel: the LEP Limit 

Before we identify benchmark slepton co-NLSP scenarios for the Tevatron, and the attendant 

limits for each, we take a moment to review the limit set by LEP2. The LEP experiments 

( 4 'We do not consider spectra with left-handed sleptons at the bottom, because the electroweak Z?-term 
splits the slepton doublet, causing the sneutrinos to be the NLSP. The collider signatures of this scenario 
have been studied elsewhere [2H [23 • 
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looked systematically for gauge mediation with slepton co-NLSPs with either prompt or 
displaced decays. The results are compiled and summarized by the LEP2 SUSY work- 
ing group [33] • Slepton co-NLSPs with prompt decays were searched for in the acoplanar 
leptons final state (OS dileptons plus missing energy). The strongest model-independent 
limit comes from direct £lr£l* r production, requiring ra^ R > 96 GeV at 95% confidence 
level. Right-handed stau production is less-constrained, > 87 GeV, because taus have 
a lower experimental acceptance, and the limit on cr depends on the neutralino spectrum, 
because of possible interference between the s and t-channels [6]. Within GGM, the limit 
from acoplanar muons sets the strongest limit on the right-handed selectron mass, because 
rrie R = rrij lR > 96 GeV follows from the flavor-independent boundary condition, discussed 
above. 

OS dilepton plus missing energy is not the most promising channel to look for slepton co- 
NLSPs at the Tevatron. This is because there are significant backgrounds due to leptonically 
decaying W~W + and ti, with cxBrsi 0.6 pb and a x Br w 0.4 pb, respectively. As we will 
review in section [4], the backgrounds are much, much lower in final states with same-sign 
dileptons or trileptons. Thus, we will focus on these more promising final states in the rest 
of the paper. The trade-off here is that they are only populated through cascade decays of 
heavier SUSY states. Thus searches in these final states will necessarily be less inclusive than 
the LEP searches, and describing their constraints will require taming the zoo of possible 
GGM spectra that contain slepton co-NLSPs at the bottom. We will turn to this problem 
in the next section. 



3 Benchmark Scenarios for Slepton co-NLSPs 

As discussed above, the Tevatron sets the strongest limits on slepton co-NLSPs through 
the same-sign dilepton and trilepton channels. These final states are produced by SUSY 
cascades ending in right-handed sleptons, but the number of leptons in the final state, and 
their kinematics, are sensitive to the details of the SUSY spectrum. This means that there 
are a very large number of possibilities in GGM, spanned by the 8-dimensional parameter 



space identified in equation 2.1 We are unable to cover this full space, but we can attempt to 



cover all the qualitative phenomenological possibilities. To that end, we have identified four 
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Figure 3: The different slepton co-NLSP spectra that form our benchmarks. 

model-independent benchmark scenarios for slepton co-NLSPs at the Tevatron. By studying 
them and optimizing experimental searches for them, one will (hopefully) ensure that no 
large qualitative gaps remain in the experimental coverage of the GGM parameter space. 

Our simplifying approach is to consider spectra with one type of production mode, which 
then decays to right-handed sleptons. The possibilities are shown in figure [3] wino produc- 
tion, higgsino production, and left handed slepton production with and without an inter- 
mediate bino. For each benchmark scenario we choose parameters, shown in table [TJ that 
maximize the branching ratio to leptons I = e, /i, which are more experimentally constrained 
than taus. We stress that we have chosen the minimal spectra that produce same-sign 
dilepton and trilepton final states. It is possible to combine our benchmarks with addi- 
tional heavier states that produce more elaborate decay chains. Our results correspond to 
conservative limits on such scenarios. 

3.1 Wino 

The first benchmark we consider consists of winos that are heavier than right-handed slep- 
tons. Since there is no coupling between charginos and right-handed sleptons, the charged 
wino always decays to fi, 

W ± -^f^ + u T (3.1) 
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(I) Wino prod. 


(II) Higgsino prod. 


(III) Slepton prod. 


(IV) Slepton prod. + 
intermediate bmo 


m~ eL = m fL 


1000 


1000 


vary 


vary 


TYlg R 


vary 


vary 


vary 


vary 


Mi 


300 


300 


300 


(rrie L + mcJ/2 


M 2 


vary 


1000 


1000 


1000 


fi 


1000 


vary (/i < 0) 


500 


500 


tan 


2 


2 


2 


2 



Table 1: Benchmark scenarios for slepton co-NLSPs with flavor democratic decays. Masses 
are specified in GeV. We choose a low-scale for SUSY breaking, y/~F = 50 TeV, giving prompt 
decays to the gravitino. 



Meanwhile, the neutral wino can decay down to all three generations of right-handed sleptons: 

W°^e% + e T , fi% + fi T , ff + r T (3.2) 

The decays to the first two generations happen only through mixing with the bino. Decays 
down to f\ can be mediated in the same way, through mixing with the higgsino, or through 
stau mixing. The latter two are proportional to the tau Yukawa coupling, so they are 
enhanced at large tan j3. By varying this and other parameters, we can make the neutral 
wino branching fractions to the right-handed sleptons flavor democratic, or tau-dominated. 
For the benchmark scenario, we choose parameters, shown in table [TJ that lead to democratic 
decays: low tan (3 and M 1 <C \i. 

We see that charged wino pair production, W + W~ , leads solely to an opposite-sign ditau 
final state, and not same-sign dileptons or trileptons. These final states do result from the 
associated production of a neutral and charged wino, Vt^°Vt^ ± . For this benchmark, 2/3 of 
Vt^ Vt^ =l= events have two OSSF leptons and t*, while 1/3 of the events have three r's. 

One simple way to modify this benchmark would be to add an intermediate neutralino 
that is mostly bino, mj < Mi < M%- Then, the mostly-wino neutralino can decay to the 
mostly-bino neutralino plus an on or off-shell Z . In most of the parameter space that will be 
of interest at the Tevatron, this decay is three-body and subdominant to the decays discussed 
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above, and we therefore choose not to consider this situation separately. 

3.2 Higgsino 

Now we consider a benchmark with higgsinos heavier than the right-handed sleptons. Recall 
that a spectrum with mostly-pure higgsinos consists of two neutral higgsinos, i?i,2, and 
a charged higgsino, H ± , all nearly degenerate with one another. Similarly to the wino 
benchmark, the charged higgsino can only decay to f\ plus a neutrino, whereas the two 
neutral winos, i?i 2 > can decay to e^e T , ji%[i T , or t^t t . The main difference between 
higgsino and wino production is that higgsinos include a new production mode of two neutral 
higgsinos, HfH®, which results in a 4-lepton/tau final state. We choose parameters, shown 
in table [TJ that yield flavor democratic decays of the neutral higgsinos: low tan (5 and low 
Mi. The production of a neutral plus charged higgsino results in the same final states as 
the wino benchmark. When two neutral higgsinos are produced, 4/9 of the events have 
two OSSF lepton pairs, 4/9 have one OSSF lepton pair and two r's, and 1/9 of the events 
have four taus. Similarly to the discussion for the wino benchmark, we could augment the 
spectrum with an intermediate bino, but the SUSY cascades will bypass the bino for most of 
the parameter space of interest for the Tevatron, yielding signals similar to this benchmark. 

3.3 Left-Handed Slepton 

The third scenario we study consists of left-handed sleptons that are produced and decay 
to right-handed sleptons. The production modes are Z^z>, uu*, and Recall that the 

left-handed sleptons and sneutrinos are split by the electroweak D-terms, with the former 
always slightly heavier, 

m~ L — m? = |cos2 (3 '{m^ (3.3) 

We will now discuss how II and v decay. The decays are quite different for the first two 
generations and for the third generation, so we will discuss them separately. 

The first two generation left-handed sleptons can three-body decay to the right-handed 
sleptons via an off-shell neutralino, 

t L -> e| + e T + r, + ^ + ft + r^ + e- (3.4) 
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Depending on the -ino masses and tan /3, these decays can either be flavor democratic or 
tau-rich. Again, for our benchmark scenario, we choose parameters (shown in table [T]) such 
that these decays are flavor democratic. We also see from (3.4) that the decay of can 
lead to either a SS or OS lepton pair. The relative branching fraction depends on the ratio 
of the off-shell neutralino mass to the slepton mass, with the SS dominating as this ratio is 
increased [32] . 

The above decays typically dominate, except when the left-handed and right-handed 
sleptons are squeezed or when the neutralinos are too heavy. In these cases, the left-handed 
sleptons prefer to decay to the sneutrino through an off-shell W, 

£ L ^D e + W~* (3.5) 

The sneutrino decays are relatively simpler. For the first two generations, the only 
possibilities are again three body decays through off-shell neutralinos: 

v t -)> e^ + e T + ui, + [i T + v e , ft+r^ + vt (3.6) 

As above, depending on the -ino masses and tan /3, these decays can either be flavor demo- 
cratic or tau-rich, and our choice of parameters in table [T] leads to flavor democratic decays. 

Finally we come to the third generation left-handed sleptons. While three-body decays 
analogous to Eq. [3T4 are also available for the left-handed stau, f 2 , the stau can also decay 
through tau-mixing and an on or off-shell Z, 

f - _> f - + Z W (3.7) 

For the parameters of table [TJ this is the dominant decay mode until f\ and fjj become 
squeezed, and then Eq. |3.5 dominates. Also, for these benchmark parameters, the third- 
generation sneutrino prefers to decay through stau mixing via an on-shell or off-shell W, 

v T ->■ rf + W+W (3.8) 
Evidently, the third generation sleptons dominantly produce OS tau pairs, as can be seen by 



inspecting equations 3.5, 3.7, and 3.8 (Additional leptons can come from the W and Z's, 
but then one must pay the branching fraction price, and plus these are generally softer.) 
Due to the lower experimental sensitivity to taus, we find that the production of the third- 
generation sleptons can be ignored when setting limits, relative to the production of the first 
two generations. 
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3.4 Left-Handed Slepton with Intermediate Bino 

The left-handed slepton decays are significantly simpler if there is a bino between the left- 
handed and right-handed sleptons. This situation deserves special consideration and con- 
stitutes our fourth and final benchmark (we pick parameters shown in table [T]). Now, the 
left-handed slepton states, II and u, all decay down to the bino plus their superpartner, and 
the bino decays flavor democratically and charge democratically to the right-handed slep- 



tons. This means that the three-body decays through off-shell Z and W (Eqs. 3.5 3.7, and 



3.8 ) are bypassed, resulting in more leptons in the final state and, as we will see, a stronger 



experimental limit. 



3.5 Benchmark Collider Signals 

We conclude this section by giving a brief characterization of the multilepton collider signals 
of the four benchmark models. Table [2] summarizes the above discussion, listing the relevant 
production modes and decay topologies for each benchmark. Here we are assuming that the 
spectra are unsqueezed, but other than this assumption, the topologies are independent of 
the NLSP and production masses. The final two columns show the percentage of events, 
weighted by the a x Br of each topology, that contain same-sign dileptons or trileptons, for 
m eR w 150 GeV, m pro d ~ 200 GeV. (m pro d denotes the mass of the particles being produced: 
wino, higgsino or slepton.) These leptons (again, e or /x) were required to have \rj\ < 1 and 
Pt > 10 GeV, and were simulated using the default PGS4 CDF simulation (so they include 
crude lepton ID cuts and parametrization of central calorimeter cracks) [28J. 

We comment that our strategy of choosing simple spectral types is distinct from the 
approach, advocated elsewhere [31], of enumerating the limits on all possible final state 
topologies (determined by the final state multiplicities and kinematics). This is because 
even the simplest possible spectra can lead to many final state topologies at once, with 
different branching fractions. Our four benchmark scenarios together contain 32 distinct 
final state topologies, and it would be cumbersome to present separate limits for each. 

Figure [4] contains histograms of the number of leptons produced for each benchmark, 
subject to the same cuts, and for the same choice of parameters. We see that for wino 
and higgsino production, the fraction of events with more than two leptons is much smaller 
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Benchmark scenario 


T~> J J 

rrod. mode 


lopologies 


7o oS Dilepton 


\ O T J. 

> 3 Leptons 


Wino 


± 

Xi Xi 


2*: It, 0t3r 


8.8 


5.3 


Higgsino 


± 
Xl Xl,2 

Xh°2 


20r, 0£3r 
4£0r, 2£2r, 0£4r 


16.5 


11.9 




vm 


4£0r, 2£2r, 0£4r 










5£0r, 3£2r, l£4r 






Slepton 


7)4- JfTL 

V T V T 
r 2 r 2 


6£0r, 4£2r, 2£4r 
0£ 2r (+2W) 
0£2t (+W + Z) 
0£2r(+2Z) 


38.7 


31.4 






4£0t, 2£2t, 0£4t 










5£0t, 3£2t, UAt 






Slepton-bino 


V T V T 
'2 '2 


be Or, 4-ezr, z^4r 
4£0r, 2£2r, 0£Ar 
UIt, 2£3r, 0£5t 
4£2t, 2£4t, 0£Qt 


51.5 


43.0 



Table 2: Production modes and final state topologies, for the flavor-democratic benchmark 
scenarios. In the final two columns, we show the fraction of events for each benchmark that 
contain same-sign dileptons or trileptons, inclusively, for m NLS p ~ 150 GeV and m prod rs 
200 GeV. Here the leptons (meaning, e or /i) are required to satisfy \rj\ < 1 and p T > 10 GeV, 
and are simulated with the default PGS4 CDF simulation. Note that the benchmarks with 
slepton production have topologies with more leptons than the benchmarks with -ino pro- 
duction, resulting in a higher branching fraction to same-sign dileptons and trileptons. 

than for slepton production. We also notice that for slepton production, events with 4 or 
5 leptons are not uncommon. (Including heavier states in the SUSY cascade, e.g. winos or 
higgsinos, would lead to even more leptons, up to 8 in all.) Finally, we learn that for slepton 
production without intermediate bino there are many more zero lepton events compared to 
slepton production with intermediate bino. This is a sign of the lepton degradation in the 
third generation production modes, as discussed above. 
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Figure 4: Number of leptons per event for each benchmark, simulated with PGS and requiring 
\r]\ <l,p T > 10 GeV. 




Figure 5: Fj t distributions for the four benchmark scenarios, with 100k events. In all 
cases, we have taken m proc i = 200 GeV. The blue (red) curves are for rug = 100 GeV 
(m^ = 150 GeV). The solid curves in the left and right plots correspond to wino and slep- 
ton production, respectively, and the dashed curves correspond to higgsino production and 
slepton production with intermediate bino, respectively. We see that all four benchmarks 
have significant J£ T > 50 GeV. 



In fig. [5j we show the missing energy distributions (again simulated using PGS) for the 
four benchmark scenarios, and two choices of the m; mass. We see that models with 
slepton co-NLSP, with masses in the range of interest, are characterized by significant J£ T > 
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50 GeV, ( 5 ) This should be compared with the typical J£ T distribution in the mSUGRA 
scenario (see appendix [C]), which is generally much softer. 

Finally, in fig. [6j we show the p^s of the three hardest leptons, for the wino and slepton 
benchmarks, again for the same choice of parameters. The higgsino and slepton with bino 
benchmarks are not shown, but the distributions look similar to the wino and slepton bench- 
marks, respectively. We have again taken m| = 150 GeV and m pro d = 200 GeV. We have 
broken down the distributions into leptons coming from Br, JIr decay (red), leptons coming 
from t decays (yellow), and leptons coming from decays of heavier states (green). We see 
that the hardest lepton tends to be quite energetic > 50 GeV, and mostly composed of 
direct right-handed slepton decays. The composition shifts towards r decays and decays of 
heavier states as one moves to the second and third hardest leptons. We note that the third 
hardest lepton is much more energetic and numerous in the slepton production scenario, 
compared with the wino scenario. 



4 Overview of Tevatron Searches 

We now review the existing Tevatron searches which are most relevant for GMSB with 
slepton co-NLSPs: same-sign dileptons+ Fj t and trileptons+Fj T . We will use these searches 
to determine the current limits and future reach on slepton co-NLSPs in section [5} 



4.1 SS Dileptons 

The simplest leptonic channel with low SM background is two leptons of the same charge. 
Same-sign dileptons have been searched for at CDF [H] and DO [13] . both with 1 fb^ 1 . The 
CDF search uses an inclusive approach and includes the ee, e/x and \x\x channels. The DO 

We also learn that for higgsino and wino production, the J5 T distribution is largely independent of 
mi . This surprising fact is because in these cases, the JE T comes mainly from the (vector) sum of the 
gravitino momenta plus the neutrino momentum coming from xt + v T ■ Since the neutrino momentum 

is ~ m proc i — rriji R while the gravitino momenta are ~ m^ R /2 each, the total J? T is ~ m proc i on average, i.e. 
the right-handed slepton mass cancels out. For the slepton production modes, there are also neutrinos from 
the decays of the left-handed sleptons, but they do not take up the entire energy, so the cancellation is less 
prominent. 
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Figure 6: The three hardest lepton p^'s for the wino and slepton benchmarks, with 100k 
events simulated with PGS and requiring \rj\ < 1. We have fixed m| = 150 GeV and 
(or mi ) = 200 GeV. The red curve represents leptons coming from cr and p,R decays. The 
yellow curve is leptons coming from r decays. The green curve is leptons coming from decays 
of heavier states in the SUSY cascade. The blue curve represents the sum of all these curves. 



search includes only the channel, and the cuts are not inclusive but are instead tuned to a 
particular mSUGRA-type model point. For this reason we choose to focus here on the CDF 
search. We also comment that CDF has conducted another same-sign dilepton search with 
2.7 fb _1 [35J. However, motivated by fourth generation quarks, they look in the / ± / ± 6jp T 
channel, in particular requiring at least two jets, one of which must be tagged as a 6-jet. 
This search is therefore not useful for constraining our benchmark scenarios and we will not 
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SS Dilepton Cuts 



Vt,Vt > 20 > 10 GeV 

|^7l,2 1 < 1 

m 12 > 25 GeV 



Table 3: The cuts employed by CDF for the 1 fb 1 same-sign dilepton search 



consider it further 



(6) 



The CDF search employs a simple set of inclusive cuts on like-sign dilepton pairs, shown 
in table [3} With these cuts, 44 events are observed in the data, with an expected background 
of 33.2 ±4.7. The largest backgrounds are Drell-Yan plus an untagged photon conversion and 
leptonic decays of dibosons, WZ and ZZ. Taus are not specifically tagged, but the search is 
sensitive to leptonically decaying taus, which account for 35% of tau decays. The strength 
of this search is its inclusive approach, which allows for a straightforward interpretation of 
its limits on GMSB. 

The CDF paper provides detailed kinematic distributions, such as E T and lepton p^'s. 
These allow us to optimize somewhat for slepton co-NLSPs. The simplest option is to 
augment the inclusive CDF cuts with a hard MET cut of E T > 60 GeV. This cut alone 
reduces the data to 4 events and the background to 2.2 events, while leaving most of the 
GMSB signal (as can be seen from fig. [5]). In the rest of the paper, we will use this selection 
to set limits and infer reach for slepton co-NLSPs. 



4.2 Trileptons 

Both Tevatron experiments have also searched for new physics in the trilepton channel, which 
has low SM background. The most recent CDF results [15] use 3.2 fb -1 of data; these are a 
straightforward update of an earlier 2 fb -1 analysis [16]. Meanwhile, the latest DO results [17] 
use 2.3 fb -1 . The CDF search uses a cut-based analysis that is simple to reproduce with 
PGS. On the other hand, the DO search uses a neural network to tag hadronic r decays, and 

( 6 ) Figure 1 of Ref. [3 5) does show the missing energy distribution of same-sign dilepton events without 
requiring jets in the event. However, the uncertainties here are 0(100%), so this figure does not appear to 
be useful for setting limits. 
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Inclusive Cuts 


Non-Inclusive Cuts (Vetoes) 


III 


Pt,Pt,Pt > 15,5,5 GeV 


Charge 


SQ = ±3 


UT 


VtiVtiVt > 15,5,10 GeV 


Z 


76 < mm- < 106 GeV 


M< i 

F, T > 20 GeV 


4th lepton 


p T > 10 GeV and \r]\ < 1 


2nd jet 


E T > 15 GeV and \r]\ < 2.5 



Table 4: The cuts employed by CDF for the 3.2 fb _1 trilepton search [15]. The inclusive 
cuts on the left select events with three leptons or two leptons and a track. Then, events 
are vetoed if the lepton charge sums to ±3, if the leptons reconstruct the Z, or if there is a 
fourth lepton or at least two jets. 

this makes it more difficult to simulate the search in order to interpret the limits on GMSB. 
For this reason, and because of the larger luminosity, we choose to focus on the CDF search 
of Ref. [If] . 

We now give a simplified description of the CDF trilepton search, focusing on the cuts 
that are important for understanding the limit on slepton co-NLSPs. For a full description 
of the analysis, we refer the reader to Refs. [16], [T7] and especially the thesis [32] • The CDF 
trilepton search selects events from two categories: three leptons, III, or two leptons and one 



isolated track 117 ^ The latter category is designed to catch 1-prong hadronic r decays, 
and this significantly improves the acceptance for models with r's in the final state, albeit 
with higher backgrounds. The cuts of the search are listed in table [4} The trileptons must 
satisfy the inclusive cuts shown to the left, and events are vetoed if they satisfy any of the 
criteria on the right. 

After cuts, one event is observed in the III category, with 1.5 ±0.2 expected; and 6 events 
are observed in the UT category with 9.4 ±1.4 expected. The main backgrounds to this 
search include leptonic decays of dibosons WZ/ZZ and ti, and Drell-Yan accompanied by 
a photon conversion, isolated track, or fake lepton. The low background rates are why the 
trilepton channel places sensitive limits on new physics with lepton-rich final states. 



( 7 ) Leptons are further broken down into tight and loose categories, where the tight leptons pass stricter 
ID requirements. The distinction between these lepton types will not be very important here, and for our 
simulation we sum the tight and loose bins, as discussed in appendix p£| 
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Although we have chosen to focus here on the CDF trilepton search, we comment that 
the DO search may be more sensitive to some GMSB models. It has no jet veto for some 
of its search bins, and it allows up to two hadronic 1-prong taus, which means it probably 
has better reach for GMSB models with colored production or with r-rich decays. We will 
revisit this issue in section |6j 

5 Results 

In this section, we present the current Tevatron limits, 95% exclusion reach, and discovery 
potential for the four benchmark slepton co-NLSP scenarios described in section |3j As 
discussed in section [4], we will focus on the CDF trilepton and same-sign dilepton searches, 
as these should provide the best sensitivity to slepton co-NLSPs. (Keep in mind, for the 
latter we are including a J£ T > 60 GeV cut as inferred from the plot of Fj T in [E].) Since these 
searches did not discuss slepton co-NLSPs, we have estimated on our own, using standard 
Monte Carlo tools, the signal acceptances for the various analyses. For more details on our 
simulations, see appendix [A} 

5.1 Current Limits and Estimated Exclusion Reach 

The current 95% limits on our four benchmark slepton co-NLSP scenarios are shown in fig. [7j 
In addition, we show the expected limit that could be set with 10 fb _1 of data. Since this 
plot is among the main results of our paper, let us take a moment and highlight some of its 
features. 

For each of the four benchmarks, we find that the limit depends primarily on the mass 
of the heavier state that is being produced (wino, higgsino, or left-handed slepton). This is 
consistent with the fact that the production cross-section is a steep function of mass, so 0(1) 
variations in the signal acceptance do not result in big changes in the exclusion contour. We 
see that for the existing limits, the constraint from trileptons is always stronger because it 
enjoys a larger luminosity (3.2 fb _1 ) than the latest dilepton search (1 fb _1 ). However, the 
projected limits at 10 fb _1 show dileptons winning out decisively. 

Of course, this is not an entirely fair comparison, since we have used the J£ T distribution 
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Figure 7: The current (filled) and projected (dashed) 95% exclusion limits, from the trilep- 
ton [15J (blue) and same-sign dilepton [Hj (red) searches, on the four benchmark scenarios 
of table [l| To calculate the limits, we have used the CL S statistic [37] , which is designed to 
limit the impact from downward fluctuations of the background. We also show the limit on 
the right-handed slepton masses (green band), as set by LEP2 [33] . 



in [13] to slightly optimize SS dileptons for slepton co-NLSPs. We believe that a similar 
improvement is possible for trileptons if the cuts are optimized (for our suggestions in this 
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channel 


m pTod [GeV] 


a[fb] 


trileptons 


bb dilepton 


111 

III 


UT 


;±;± 
FP 


wino 


190 


69 


8 x 10~ 3 


.02 


.04 


higgsino 
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47 


.02 


.02 


.08 


slepton 


155 


41 


.02 


.01 


.11 


slepton +bino 


170 


23 


.04 


.03 


.24 



Table 5: The masses, cross-sections, and acceptances at the edge of the current 95% limits 
for the four benchmark scenarios. The cross-sections are leading order, and in each case we 
fix m~ eR = 120 GeV. 

direction, see section |7|). To emphasize this point, we note that the difference in reach 
between the trilepton and same-sign dilepton searches is largest for the two left-handed 
slepton production benchmarks. This is because left-handed sleptons lead to higher lepton 
multiplicities, where the same-sign dilepton search benefits from its inclusive approach, and 
where the trilepton search loses out by vetoing on additional central leptons (among other 
things) . 

It might seem surprising that the projected 10 fb _1 limit on left-handed sleptons is so 
strong, since their cross-section is a factor of 5-10 smaller than the wino or higgsino produc- 
tion rates, as shown in figure [2] But again, this is because left-handed slepton production 
leads to more leptons in the final state (see table [2] and figure [4]), and therefore a higher 
experimental acceptance from SS dileptons. The signal acceptances of the various searches, 
calculated at the boundaries of the excluded regions, are summarized in table [5j 

The experimental sensitivity of the trilepton and SS dilepton searches can degrade when 
the mass splitting between the production mode and right-handed sleptons is small. The 
reasons for this vary between the different searches and the different benchmark scenarios. 
For winos, the trilepton search has a squeezed region, simply because the third lepton p T 
is being squeezed out. The dilepton search has no analogous squeezed region, since there 
the two SS leptons can come from the decay of the co-NLSPs. For higgsino production, the 
trilepton and dilepton searches have no squeezed region. This is because there is always a 
sizeable mass splitting between the two neutral higgsinos, — ~ 10 GeV, so the 
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third lepton can come from the decay of mg 2 to the co-NLSPs. Finally, for the two slepton 
production scenarios, both the trilepton or dilepton searches have a squeezed region. Here 
the reason is because the sneutrinos decay to right-handed sleptons through a three-body 
decay in the squeezed regime, but for low-scale gauge mediation, this three-body decay can 
lose to direct sneutrino decay to gravitino, v — > v + G. Therefore, in the squeezed regime, 
the sneutrinos become co-NLSPs and there is not enough branching ratio into trilepton or 
like-sign dilepton final states to produce a limit. 

5.2 Estimated Discovery Potential 

Finally, let us also estimate the potential for 3cr and 5a discovery in 10 fb -1 . Here we spe- 
cialize to a model line, so that we can show the confidence level statistic in more detail. 
The result is shown in figure [SJ fixing = 120 GeV. For each channel, we find a signifi- 
cant range of parameters where 3a evidence is possible with the same-sign dilepton search. 
Specifically, 3a evidence is possible up to 215 GeV for winos, 235 GeV for higgsinos, and 
230 (220) GeV for sleptons with (without) an intermediate bino. There is even a signifi- 
cant range of left-handed slepton masses consistent with the current limit (about 30 GeV 
in mass), where a 5a discovery is possible at the Tevatron! It is exciting to see that such 
a dramatic discovery is still possible at the Tevatron. This significant discovery potential 
is possible because the current searches have not been optimized for the slepton co-NLSP 
type signature, and because the same-sign dilepton search has not been updated since the 
relatively small luminosity of 1 fb -1 . 

6 Beyond Our Benchmarks 
6.1 Tau Rich Decays 

Above, we have explored the Tevatron limits and reach for gauge mediation models with 
flavor democratic decays. This is because the existing Tevatron searches for multilepton 
events have the best sensitivity in this case. But as discussed in section |2j it is also a 
very common situation to have decays producing mostly rs. This can happen either with 
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Figure 8: The 10 fb^ 1 discovery potential of the trilepton [15] (blue) and same-sign dilep- 
ton [Hj (red) searches, fixing mg H = 120 GeV. We use the CL;, statistic [37J which measures 
the probability of the background fluctuating to give the signal plus background. The green 
(purple) lines show the thresholds for 5a (3a) discovery and the blue (red) shaded region 
is already excluded at 95% by the existing trilepton (dilepton) search. For all four bench- 
marks we find that both searches can discover supersymmetry at 3a for a range of allowed 
parameters, and the SS dilepton search has the potential for 5a discovery for all but wino 
production. 



slepton co-NLSPs but preferential decays to f\\ or with stau NLSPs, where all SUSY cascades 
necessarily end in f\. 

Such models are obviously much less constrained at the Tevatron than models with 
slepton co-NLSP and flavor democratic decays. The existing Tevatron searches are mainly 
sensitive to these events when the rs decay leptonically (which occurs with a 35% branching 
fraction). Since the leptons arise from a three-body decay r — > t + vg, + u T , the daughter 
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Figure 9: The current limit (solid) and future reach (dashed) from the CDF same-sign 
dilepton and trilepton searches, for wino production with varying branching fraction to fi 
and wino mass. For this plot, the right-handed slepton masses were fixed at the LEP II 
limit, namely 96 GeV. The right-hand side of the plot corresponds to the flavor democratic 
case studied above. Moving leftward on the plot increases the fraction of events populated 
with rs, worsening the experimental sensitivity. 

lepton is significantly softer than the mother r. The poor experimental acceptance for final 
states with r's arises from a combination of branching ratios and softer leptons. Note that 
the trilepton search will also be sensitive to events with one-prong hadronic r (50% branching 
fraction), provided the other r's decay leptonically. Here however the backgrounds are much 
higher. 

To illustrate how the sensitivity of the Tevatron searches degrades when models are tau 
rich, we show the 95% limit and reach contours on figure |9j for a simple deformation of 
the wino benchmark scenario defined in section [2j Recall that in the original benchmark, 
tan/3 = 2 and Mi = 300 GeV, and so the neutral wino decays democratically to fi, en, and 
£lr. In fig. |9j we have considered a modified scenario where Br(W/° — >■ f\ + r) is allowed to 
vary. (One can think of this as varying tan /3 and Mi.) As can be seen from the figure, the 
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reach and limit are significantly degraded as Br(W /0 — > fi+r) — > 1. In fact, there is currently 
no limit whatsoever from Tevatron searches for the fully r-rich case. The trilepton search 
loses sensitivity at lower values of because the third lepton, which must come from 

wino decays, becomes too soft. However, at higher values of the wino mass, trileptons have 
the best reach, because they include the UT category This should be contrasted with the 
flavor democratic case, where SS dileptons have the best reach, because of its more inclusive 
nature. 

The lack of any current Tevatron limit on tau-rich scenarios is a huge gap in the search 
for GMSB. This may present an opportunity for early LHC running if dedicated searches are 
developed for tau-rich final states. More generally, the results presented in this subsection 
highlight the importance of developing search strategies that are more sensitive to hadronic 
taus. We conclude this section with an important caveat to the above discussion. We have 
not studied the limits resulting from the DO trilepton search [17J, which was conducted with 
2.3 fb _1 and includes a search category with a muon and two (hadronic 1-prong) taus. This 
bin may improve on the limit presented in figure [9j One difficulty for deriving such a limit 
is that the DO search relies on a neural network for tagging hadronic tau decays, which is 
difficult to replicate using PGS. 

6.2 A Look Towards the LHC 

In this paper, we have focused on the limits and reach, resulting from electroweak production 
at the Tevatron. The next step will be to search for gauge mediation with slepton co-NLSPs 
at the LHC (see for example the multilepton signatures discussed in [2S1 ESI ED]), which is 
now running with center of mass energy of 7 TeV. In order to understand how well the LHC 
can compete with the Tevatron, we can compare the ratios of cross-sections of interest. For 
example, 200 GeV left-handed sleptons satisfy o~lhc I 'o~T e v — O(10). More generally, the 
Tevatron with ~ 10 fb _1 will produce at least as many light uncolored states as the LHC 
with < 1 fb _1 . Meanwhile, many of the backgrounds will be larger at the LHC, implying 
that the Tevatron retains an advantage for electroweak production |18j . On the other hand, 
the early LHC enjoys a significant advantage for gg initiated processes, and 450 GeV gluinos 
satisfy o~lhc I o~TeV — 500, where we have fixed the squark masses to 1 TeV. Colored SUSY 
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production at the LHC will very quickly overtake the Tevatron, providing natural channels 
to look for GMSB in the first run of the LHC. We stress that the early LHC and the 
Tevatron are largely complementary experiments from the point of view of SUSY discovery - 
the former is well-suited for light colored states, while the latter (owing to its large integrated 
luminosity and lower QCD backgrounds) is ideal for light EW states. 

We will address colored production of slepton co-NLSPs at the early LHC in an upcoming 
publication |30j. As a preliminary step, we here establish the limit from the Tevatron for 
an example process. We consider slepton co-NLSPs with a low energy spectrum consisting 
of a light gluino, g, flavor-degenerate right-handed sleptons, Ir, and an intermediate bino, 
trig > m B > m^S 8 ^ The gluino decays to the bino through a three-body decay mediated by 
an off-shell squark, 

g^B + 2j (6.1) 

and then the bino decays democratically to the right handed sleptons, as in our fourth 
benchmark above, 

B^(ei,e*), (/4/i^), (f±r T ) (6-2) 

with equal branching ratios (1/6 per final state). And as usual, the slepton co-NLSPs decay 
to their superpartner and a gravitino. The resulting final state is 4/ + 2j + fir 



(0) 



For 



this scenario, we choose a low SUSY breaking scale with a gravitino mass of tuq = 1 eV, 
and we fix the bino mass to be halfway between the gluino and right-handed sleptons, 
m B = (m g + mj )/2. When determining the gluino production cross-section, we assume 
degenerate squarks of mass irtq = 1 TeV, and we use the NLO value of Prospino [4"T1 02] . 
We present the current limits from the CDF same-sign dilepton and trilepton searches, 



along with the 10 fb" 1 reach, in figure 10 In table [6j we show the gluino mass, production 
cross-section, and experimental acceptances at the edge of the current 95% exclusion region. 
We see that the dilepton search has an acceptance that is 10 3 — 10 4 times larger than the 
trilepton search. This is because the CDF trilepton search includes a jet veto that throws 
out most of the signal (this was designed to reduce the background from tt). Consequently, 



( 8 'Keep in mind that such a spectrum is perfectly allowed in GGM, but never occurs in MGM where the 

gluino is always much heavier than the bino and the right-handed sleptons. 
( 9 ^A similar final state is produced by a spectrum with the bino lifted, rrig > rrig, in which case the gluino's 

dominant decay is 4-body, but we will not consider this situation here. 
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mg [GeV] 


a [fb] 


trileptons 


SS dilepton 


III 


UT 




400 


36 


10" 4 


4 x 10" 5 


.19 



Table 6: The gluino mass, production cross-section, and acceptance at the edge of the 
current 95% limit. We fix m e - fl = 120 GeV, and the cross-section is the NLO value from 
Prospino [HI EE]. 
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Figure 10: The current limit (solid) and future reach (dashed) for gluino production with 
slepton co-NLSPs. We use the CDF same-sign dilepton and trilepton searches, and choose 
a benchmark scenario with an intermediate bino with mass halfway between the gluino and 
right-handed sleptons. The trilepton search includes a jet veto and therefore produces a 
weaker limit than the same-sign dilepton search, which is inclusive. 

we find that the same-sign dilepton search sets a significantly stronger limit than the CDF 
trilepton search, despite the smaller luminosity of 1 fb -1 versus 3.2 fb -1 . We comment that 
the DO trilepton search [17] did not use a jet veto for some of their search bins, and we 
therefore expect this search to set a stronger limit that the CDF trilepton search, for this 
scenario. We do not derive this limit here. It would also be interesting to infer a limit from 
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the CDF search for / ± / ± 6j^ T [35J, although their requirement for a jet with a 6-tag makes 
them sensitive only to SUSY cascades involving heavy flavor squarks. 

7 Discussion 

7.1 How to Optimize for GMSB 

In this section, we discuss how the above searches can be optimized for the slepton co-NLSP 
signal, and we suggest a few additional searches that could improve the Tevatron's sensitivity 
for GMSB. Our suggestions here are mostly culled from the preceding analysis. Some are 
rather obvious, but we hope it will be useful to collect them all here in one place. 

Many of the SUSY searches at the Tevatron, such as the trilepton search described 
above, are specifically optimized for the mSUGRA model. In order to understand how 
the cuts can be improved for the slepton co-NLSP scenario, it is helpful to compare the 
multilepton signals of mSUGRA and GMSB. The details of this comparison are contained 
in appendix [Cj In summary, we learn that slepton co-NLSPs produce more leptons, more 
missing energy, and harder leading lepton pr's. Therefore, our main point is that when 
optimizing for gauge mediation one should choose cuts that are as inclusive as possible 
regarding the lepton multiplicity. Backgrounds can be reduced while keeping the signal by 
cutting hard on missing energy and lepton pr's. 

1. Same-sign dilepton optimization: The same-sign dilepton search can be optimized 
for GMSB by taking advantage of the hard MET and hard lepton px spectra of the 
slepton co-NLSP scenario. We find that the search gives significantly stronger limits if 
a hard MET cut is added, E T > 60 GeV, and we used this approach when presenting 
the limits and reach in section [5j Moreover, it may be possible to further optimize the 
search by cutting tighter on lepton pr's. 

2. Trilepton optimization: The current CDF trilepton search loses significant accep- 
tance for slepton co-NLSPs because of the various vetoes it includes. For example, the 
fourth lepton veto and requirement that the trilepton charge not sum to ±3 removes 
many GMSB events with more than 3 leptons in the final state. In addition, the Z 
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veto removes a significant fraction of signal events where pairs of leptons accidentally 
reconstruct the Z. The cumulative effect of these vetoes degrades the trilepton search 
acceptance, which is the main reason that we found significantly better reach with the 
more inclusive approach of the same-sign dilepton search, in section [5j The trilepton 
search also includes a jet veto, which we found to reduce the acceptance for slepton 



co-NLSPs coming from colored SUSY production (see figure 10 and table [6] above). 
Our suggestion is to use a more inclusive approach, by removing the lepton, jet, and Z 
vetoes and the charge requirement. Instead, backgrounds can be reduced by requiring 
additional MET and harder first and second lepton p^'s. 

3. > 4 lepton search: To our knowledge, the Tevatron has not searched for new physics 
in the > 4/ + E T channels, although such searches were suggested by some of the early 
GMSB literature IH0IE]. ft would be very useful to conduct searches in these channels 
because they have low SM background and GMSB can produce many more than three 
leptons (as many as 6 with our benchmarks, see tableland figure El) in the final state. 
Such a dedicated search would be particularly effective for our slepton benchmarks, 
where the production of > 4 leptons is common. 

4. Multi-tau search: As we found in section [6j the same-sign dilepton and trilepton 
searches have poor acceptance for slepton co-NLSP scenarios with r-rich decays. This 
is because these searches mostly rely on leptonically decaying rs, which requires paying 
a 35% branching ratio per r. Therefore, we recommend pursuing multiple hadronically 
decaying rs in the final state. While this is more difficult experimentally, it may be 
possible to take advantage of the large E T and large r pr's of the slepton co-NLSP 
scenario. 



7.2 Excess in SS Dileptons? 

The CDF same-sign dilepton search [13] , described in section |4j has looked inclusively for 
events in 1 fb _1 with two leptons of the same charge. Interestingly, the inclusive search has 
a slight excess, which is barely consistent within 2a. Even more interestingly, CDF shows 
several distributions for events that pass the inclusive selection, including histograms of the 
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MET and hardest lepton p?. CDF has observed four events with MET > 80 GeV when only 
~ 0.9 are expected and four events with leading lepton p? above 90 GeV where ~ 1.5 events 
are expected. It is also intriguing that Sleuth/ Vista have identified like-sign dilepton final 
states to be among the most discrepant states, relative to background expectations, at both 
CDF [43J and DO [H]. 

These same-sign dilepton excesses are not statistically significant. Nevertheless, it is 
amusing to entertain the idea that CDF has detected a hint of supersymmetry, because, as 
we have seen in section|3j same-sign dilepton events with high MET and a hard leading lepton 
are generically predicted by gauge mediation with slepton co-NLSP. So in this subsection, 
we will attempt to answer the following question: can gauge mediation explain this slight 
excess at high MET and lepton pt, consistently with the current experimental constraints? 
(For a different possible explanation, see [4"5]). 

In order to answer this question, we have looked at models at the edge of the current 



95% exclusion region, as determined in section [5] and 6.2 and counted the number of events 
that pass the CDF same-sign dilepton cuts with high MET or hard leading lepton p?. For 
simplicity, we have fixed mg H = 120 GeV. The result is summarized in table [7j We find that 
all of our benchmark scenarios can give several events of interest. For example, left-handed 
slepton production with (without) an intermediate bino gives about 4 (3) events with MET 
above 80 GeV. We find this exciting because the trilepton searches could have already ruled 
out any potential explanation of the "excess" in terms of slepton co-NLSPs. This is not the 
case. We have also considered gluino production with slepton co-NLSP, which we discussed 



above in section 6J2 For this channel, we find that gauge mediation could have produced as 
many as ~ 6 events with high MET or a hard leading lepton, consistent with current limits. 

In more detail, we have investigated the range of currently allowed production masses 
which could "explain" the slight excess (by which we mean, rather arbitrarily, to have > 2 SS 
dilepton events with E T > 80 GeV) . For higgsino production, this corresponds to production 
masses in the range mg — 190 — 210 GeV. For left handed slepton production, we find 
ma = 150 — 190 GeV without an intermediate bino, and = 170 — 200 GeV with an 
intermediate bino. For gluino production, the mass range of interest is m~ g = 400 — 430 GeV. 

We are encouraged to learn that all of these mass ranges are within reach at the Tevatron 



with 10 fb 1 (see figures [7] and 10). In fact, these ranges of masses can be discovered with 5<r 
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channel 



wmo 
higgsino 
slepton^ 
slepton^ +bino 
gluino 



Number of Events (1 fb 1 ) 



MET > 80 GeV 



1.8 
2.8 
3.1 
3.9 
5.6 



p\ > 90 GeV 



0.9 
2.2 
2.4 
2.9 
6.8 



Table 7: The number of same-sign dilepton events with high MET or high leading lepton p?, 
for models at the edge of the current 95% exclusion (see section [5|, fixing m^. R = 120 GeV. 
Such events may explain the slight excesses observed by CDF. 

significance, as can be seen in figure [8] This means that if CDF has seen the first evidence of 
supersymmetry in the same-sign dilepton search, then the Tevatron still has the opportunity 
to make a spectacular discovery. Alternatively, this hypothesis can be readily ruled out by 
upcoming data. 
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A PGS Calibration 

For event simulation in this paper, we use Pythia6.4 [27], including showering, hadroniza- 
tion, and multiple interactions. The events are then passed through PGS4 [28J for detector 
simulation. We use Suspect [29] to calculate the superpartner spectrum and Pythia decay 
tables when available. Pythia does not include the relevant two-body decays of the lightest 
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Process 






ZW 


0.094 


0.080 


ZW (fi T > 15 GeV) 


0.086 


0.070 


zz 


0.15 


0.18 



Table 8: Background ZW and ZZ acceptances for the CDF same-sign dilepton search [T4] . 
We compare the acceptances from our PGS simulation to the published acceptances and find 
good agreement. The acceptances are normalized to leptonic ZW, / = e, /x, r. 

neutralino and three-body decays of left-handed sleptons via off-shell W, Z and neutralino. 
We use a private code for these decay tables. For the left-handed slepton decays, we have 
implemented the three-body decay results from Refs. (32j SH] . These results were originally 
for three-body squark decays, and three body decays of — > fx + e + r, but they are 
straightforward to adapt to all the cases of interest. 

In this appendix, we describe the consistency checks and calibrations of our event sim- 
ulations. We have confirmed the accuracy of our simulation of the trilepton and like-sign 
dilepton searches by comparing PGS output to the published acceptances for signal and 
some backgrounds. We find that the dilepton search is well-simulated by PGS out-of-the- 
box, whereas minor modifications are necessary to reproduce the results of the trilepton 
search. 

A.l SS Dileptons 

For the CDF dilepton search [TJ], described in section [5j we use PGS "out-of-the-box" (i.e. 
without any further modifications), and we compare the acceptance of our simulation to the 
published acceptance for the leptonic ZW and ZZ backgrounds. These are summarized in 
table[8j The cuts that have been applied are for the "inclusive" analysis described in section[5j 
For ZW where the Fj T is real, we have also compared acceptances including an additional 
Fj t > 15 GeV cut. In all cases, we see that the acceptances agree with the published values 
to better than 25%, indicating good agreement between our PGS simulations and the official 
CDF simulations. 
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e e+ 




p t > 20 GeV 


p T < 20 GeV 


p t > 20 GeV 


p T < 20 GeV 


weight 


0.80 


0.83 


0.97 


0.56 



Table 9: We define lepton weights so that PGS matches the CDF simulation in the dilep- 
ton control region. First, the weights for pt > 20 GeV are determined by dileptons that 
reconstruct the Z, and then the weights for pt < 20 GeV are determined by the remaining 
Drell-Yan events, which are dominated by lepton pairs with low invariant mass. 

A. 2 Trileptons 

It was more challenging to reproduce the CDF trilepton search [15j US] using PGS. In the 
end, we found that PGS out-of-the-box was not sufficient. It is easy to understand the 
reason for this: since the search requires three leptons, any difference between the per-lepton 
PGS reconstruction efficiency, e, and the actual efficiency of CDF, is amplified because the 
acceptance scales as e 3 . Also, the trilepton search probes softer leptons than the dilepton 
search, p T > 5 GeV instead of pt > 10 GeV. Softer leptons have a lower, more pr-dependent 
reconstruction efficiency which is not modeled correctly by PGS. We find that PGS, out-of- 
the-box, overestimates the trilepton acceptance by about a factor of 2. 

To model more accurately the CDF detector, we have modified PGS in three ways: 
(1) we use a crack parameterization that better matches the CDF geometry [19J, (2) we 
have modified the lepton definitions to closer represent the ones used by the CDF trilepton 
search [36], and (3) we have defined an isolated-track object class, as used in the CDF 
analysis to identify 1-prong r decays, 



(10) 



After including these modifications, the agreement with CDF is better, but still not good 
enough. To understand better the per-lepton acceptances, we have examined using PGS the 
dilepton control sample used in the CDF trilepton search [36]. This is populated mainly by 
Drell-Yan production, and its acceptance should be dominated by the per-lepton acceptance. 



( 10 )The CDF analysis includes two types of electrons and two types of muons, referred to separately as tight 
and loose to represent more and less stringent ID requirements, respectively |36j . We find that PGS leptons 
are not refined enough to make this differentiation. Therefore, we use one lepton ID class, which is meant 
to simulate the sum of tight and loose leptons used in the CDF analysis. 
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Nominal Point 


61 


Zt + t 


PGS 


8.5 


11.8 


PGS (reweighted) 


4.5 


8.1 


CDF 


4.6 


6.8 



Table 10: The PGS, PGS with lepton weights, and CDF trilepton counts for the nominal 
mSUGRA point defined in [16] . Here, m = 60 GeV, m 1/2 = 190 GeV, tan/3 = 3, A = 0, 
and jj, > 0. The two bins correspond to three leptons, and two leptons and one track. 




mo [GeV] 

Figure 11: The left plot shows the PGS limit and the right plot shows the actual CDF 
limit [16]. We have included the PGS lepton weights of table [9] We include a uniform 
K-factor of 1.35 to adjust the LO pythia cross-section for N 2 + C\ production, to the NLO 
cross-section. We have checked using Prospino [H] that the K-factor varies by less than 5% 
over the parameter space. 

We find that PGS overestimates the lepton acceptance relative to the CDF search. The 
difference between the PGS and true efficiency depends on the lepton px, particularly for 
low pt muons. Therefore, we define per-lepton weights for PGS, which we use to reweight 
PGS events to match the dilepton control sample. (A similar procedure is used by the CDF 
collaboration when tuning their detector simulator to control samples [36J . ) We allow the 
lepton weights to depend on lepton px by using two pt bins separated by px = 20 GeV. The 
lepton weights we use are shown in table |9j 



38 



With our modified and reweighted PGS in hand, we check the validity of our simulation by 
comparing to the published signal acceptances for the mSUGRA model studied by CDF [16J. 
The trilepton search reports efficiencies for a "nominal point" in mSUGRA parameter space, 



defined in |16j . In table 10, we compare the PGS event count to the published value, before 
and after applying lepton weights. We find that our simulation agrees to better than 20% 
after applying lepton weights. Since the weights were calibrated to the dilepton control 
sample, this provides a nontrivial check of our simulation. 

Finally, in figure [TTJ we compare the 95% limit determined with acceptances from our 
simulation, to the published CDF limit at 2 fb _1 , for a plane of mSUGRA parameter space. 
We find that our simulation reproduces fairly well the shape of the CDF limit. We comment 
that the PGS acceptance does still differ from the published acceptance by up to a factor of 
about 2 in regions of parameter space with at least one very soft lepton, pt — 5 GeV. 



B Application to Minimal Gauge Mediation 

So far, we have studied slepton co-NLSPs within the model-independent context of GGM. It 
is instructive to see how Minimal Gauge Mediation (MGM) fits within our framework, since 
MGM served as the basis for the original phenomenological literature on gauge mediation [3j 
HI El EJ [8j [9] . MGM is characterized by three continuous parameters: the SUSY breaking 
scale A, the messenger mass M, and tan /3; and two discrete parameters: the messenger 
number N and sign(/x). This is clearly a much smaller parameter space than that of GGM, 
and it greatly restricts the phenomenological possibilities. 

For instance, in MGM there are only three possible NLSPs: slepton co-NLSP, stau NLSP, 



or bino-like neutralino NLSP Slepton co-NLSPs with prompt decays are generic when 
iV > 3, tan/3 < 10, and M < 10 7 GeV p2 |32]. Larger tan/3 typically leads to stau NLSP, 
smaller messenger number leads to bino NLSP, and larger messenger masses imply a high 
SUSY breaking scale and lead to displaced decays. 



In figure 12, we show a typical MGM spectrum within the prompt slepton co-NLSP 



( n ) There are some special values of the MGM parameters with sneutrino NLSP (see for example the second 
footnote of [6]). These models are within the kinematic reach of LEP2 so are likely excluded, although this 
may not have been checked. 
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Figure 12: A typical MGM spectrum with slepton co-NLSPs, as a function of the messenger 
SUSY breaking scale, A. We fix the messenger number and masses, N = 4 and M = 100 TeV, 
and we choose tan (3 = 2 and sign(/i) < 0. We show the states that are relevant for the 
Tevatron. Trileptons and opposite-sign dileptons result from wino production, and left- 
handed slepton production. We have used SOFTSUSY jl?] to calculate the spectrum. 

regime. The masses in MGM are proportional to A, shown on the x-axis. Here we have fixed 
N = 4, M = 100 TeV, tan/3 = 2, and /i < 0. We find that the spectrum is similar to our 
type IV benchmark of section |2j with light left-handed sleptons and an intermediate bino- 
like neutralino with mass between the right and left-handed sleptons. Moreover, the lightest 
chargino and second lightest neutralino are both mostly wino-like, and are light enough to 
play a role in the Tevatron phenomenology. The colored states in MGM are too heavy to be 
produced at the Tevatron, and are not shown. 

MGM with slepton co-NLSP can populate like-sign dilepton and trilepton final states 
when left-handed sleptons are produced, as in the type IV benchmark, or when winos are 
produced. Winos will decay to the left-handed sleptons, resulting in final states that are 
similar to left-handed slepton production, augmented by additional leptons from the wino 



decay. We show the current and projected (10 fb _1 ) limits on MGM in figure 13, fixing 
M = 100 TeV, tan /3 = 2, \x < 0, and varying A and N. For N = 5, the current Tevatron 
limit is worse than the LEP II limit on direct p>Rjl R production, of mj R > 96 GeV. This 
follows from the fixed mass ratios between Ir and the heavier states, a restriction that is not 
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Figure 13: The CL S statistics for the CDF trilepton search (blue) and SS dilepton search 
(red) are shown for MGM lines with messenger numbers N = 3, 4, 5. The solid (dashed) 
lines indicate current (projected 10 fb -1 ) searches. The dotted black lines indicate the 95% 
confidence level. Included also in the plots are the LEP II limits on direct right-handed 
slepton production, which are stronger than the current Tevatron limits except for N = 3. 

present in GGM. For N = 4, LEP and Tevatron are competitive, and for N = 3, Tevatron 
sets a slightly better limit than LEP. Tevatron gains a competitive advantage at smaller 
iV because the winos are lighter, relative to the right-handed sleptons, leading to a larger 
cross-section at the Tevatron. For all three values of N, we find that the projected reach at 
Tevatron will surpass the LEP II limit. 

It is possible to set a conservative limit on MGM by considering left-handed slepton 
production alone. For this purpose, we can directly apply our results from section [5] We 



illustrate this procedure in figure [14J where we show the MGM mass line, superimposed on 
the type IV limit plot of section |5j This mass line, mj L ~ 2m ifl , applies robustly to MGM 
with tan (3 = 2, regardless of messenger number or messenger mass in the prompt regime, 
M < 10 GeV. The limit and reach on left-handed slepton production correspond to when 
the MGM mass line crosses the relevant limit and reach curves on this plot. The actual limit 
and reach, including wino production is about 10% stronger in mass, and shown on the plot. 
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Figure 14: The MGM mass line superimposed on the exclusion and reach contours for left- 
handed slepton production with an intermediate bino, reproduced from figure [7} MGM 
imposes a relation between the left-handed and right-handed slepton masses, corresponding 
to the brown line for tan /3 = 2. MGM also includes a light wino, and wino production gives 
additional events. This leads to a slightly stronger limit (orange circle) and reach (purple 
circle), as shown here for messenger number N = 4. 

C GGM vs. mSUGRA 

We include a brief comparison of the multilepton signals from GMSB with slepton co-NLSP 
to those from mSUGRA (for a review with references see pQ). This comparison is useful 
because many Tevatron searches for SUSY, such as the trilepton search described above, 
were designed with the mSUGRA signal in mind. This has led the searches to be biased in 
various ways. 



Shown in fig. 15 are the number of leptons, the Fj t , and the hardest lepton distri- 
butions, for a nominal mSUGRA point with tan/3 = 3, fi > 0, A = 0, itlq = 100 GeV and 
m i/2 — 265 GeV. This corresponds to a wino mass of ~ 200 GeV, a bino mass of ~ 100 GeV, 
and a right-handed slepton mass of ~ 150 GeV. The left-handed sleptons are heavier than 
the wino and are effectively decoupled. The only production mode for multiple leptons here 
is xtx2- These states decay with nearly 100% branching fraction (and nearly flavor demo- 
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Nlep MET (GeV) lepton p Tl (GeV) 



Figure 15: mSUGRA comparison. For mSUGRA we choose a nominal point, tan = 3, 
\i > 0, A = 0, mo = 100 GeV and mi/2 = 265 GeV. On the left we show the number of 
leptons with \r/\ < 1 and px > 10 GeV. Comparing with figure[4j we see that GMSB tends to 
produce more leptons than mSUGRA. The center and right plots compare the distributions 
of missing energy and hardest lepton px (with 100k events) between mSUGRA (blue) and 
GGM with wino production (red). For GGM we match the wino and right-handed slepton 
masses with mSUGRA, ~ 200 GeV and ~ 150 GeV. 



cratically) to the right-handed sleptons, which then decay to the bino LSP. So this scenario 
is most analogous to our wino production scenario, with the only difference being that the 
LSP is 100 GeV here instead of the massless gravitino. 



In fig. 15, we have also included for comparison the Fj t and pxi distributions for the 
wino production GMSB scenario with the same wino and right-handed slepton mass. We 
see that the lepton distribution is quite similar to the that of the wino production scenario 
(the left-most plot in fig. [4]). Of the slepton co-NLSP scenarios we have considered, wino 
production leads to the fewest number of leptons, which means that GMSB can produce 



significantly more leptons in the final state than mSUGRA. We can also see from figure 15 



that the Fj T and pxi distributions for the GMSB scenario are relatively harder, since the 
LSP is massless. This suggests that when optimizing for GMSB signatures, as opposed to 
mSUGRA signatures, one should prepare for the possibility of more leptons in the final state 
and one can cut harder on Fj t and lepton px, perhaps relaxing other cuts or vetoes as a 
consequence. For our suggestions in this direction, see section [7j 
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